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Abstract 

We consider three Ginibre ensembles (real, complex and quaternion-real) with a deformed measure 
and relate them to the known integrable systems via presenting partition functions of these ensembles 
in form of fermionic expectation values. We also introduce double deformed Dyson- Wigner ensembles 
and compare their fermionic representations with these of Ginibre ensembles. 
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Ginibre ensembles play an important role in many statistical problems. They were introduced in 1965 by 
Ginibre as non Hermitian analogues of famous Wigner-Dyson ensembles, namely, ensembles of random 
real symmetric, Hermitian and quaternion self-dual random matrices, known also as orthogonal, unitary 

■ and symplectic Wigner-Dyson ensembles. The non-Hermitian analogues are called real Ginibre, complex 

■ Ginibre and quaternion-real Ginibre ensemble respectively. 
Let us recall that both Wigner-Dyson ensembles and Ginibre counterparts are Gauss ones. This 

was enough for the problems considered by physisists working in quantum chaos. It is widely known 
that the deformation of the Gauss measure of each Wigner-Dyson ensemble makes them tau-functions of 
integrable hierarchies where deformation parameters treated as the so-called higher times. First this link 
" was established for the unitary ensemble which was identified with the ID Toda lattice tau function [12] . 

£SJ . and later also for the orthogonal and symplectic ones which were related to the so-called Pfaff lattice [7J . 

The origin for the deformation of WD ensembles were applications of these models to certain problems 
in physics (string theory and 2D quantum gravity and other application of summation over poly-angulated 
surfaces) and in mathematics (few enumeration problems). At present there are no special motivation 
to consider deformations of the Ginibre ensembles, however one may believe in future applications. 

Wigner and Dyson had physical reasons to consider orthogonal, unitary and symplectic ensembles. 
Physical systems are described by Hamiltonian operators which are Hermitian. Ginibre had no physical 
reasons and introduced non-Hermitian analogues from " mathematical" point of view, having a hope that 
they may be of use in future. Indeed at present Ginibre ensembles are in a focus of attention in many fields 
(such as quantum chromodynamics, dissipative quantum maps, scattering in chaotic quantum systems, 
growth processes, fractional quantum-Hall effect, Coulomb plasma, stability of complex biological and 
neural networks, directed quantum chaos in randomly pinned superconducting vortices, delayed time 
series in financial markets, random operations in quantum information theory, see a review article [T] for 
details.) 

I am going to show that the partition function of the real and quaternion-real Ginibre ensembles 
are related to integrable systems similar to many other matrix models [T2J, [T3], [2T], [23], [7], [5] which 
may be referred as deformed Dyson- Wigner ensembles. Here we use the so-called "large" BKP hierachy 
introduced in 4 (named by authors "charged" or "fermionic" BKP) and the "large" 2-BKP hierarchy 
introduced in [3J. The last hierarchy is rather similar to 2-DKP (Pfaff-DKP) one introduces in 200 We 
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1 Large BKP hierarchy includes the famous KP one as a particular reduction, while large 2-BKP includes the Toda 
lattice hierarchy |19| . In the present text we will refer these hierarchies as BKP and 2-BKP ones instead of the "large 
BKP" and "large 2-BKP" . Make difference between these "large" BKP hierarchy containing KP and the "small" (referred 
in [4] as "neutral") one introduced in 28 and which is a subhierarchy of the KP hierarchy. 
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shall see that real and quaternion-real Ginibre ensemble both deformed in a appropriate way may be 
related to these hierarchies where the higher times, namely, an integer L and a pair of semi-infinite sets 
of numbers t = (ti,t^, ■ . .) and s = (si,S2, • • •) plays the role of deformation parameters. We say that 
this is "double" deformed ensembles because two sets of parameters is used. 

To be more precise it is only the case where L > and s = may be related to what is called Ginibre 
ensembles. To consider more general deformations we need to introduce restricted (real, quaternion- 
real, complex) Ginibre ensembles where the space of random matrices consists of only invertable (real, 
quaternion-real, complex) matrices. Then we prove that partition functions of deformed (real and 
quaternion-real) Ginibre ensembles are BKP tau functions, where t and L > are deformation pa- 
rameters and BKP higher times. While partition functions of double deformed restricted (real and 
quaternion-real) ensembles are 2-BKP tau functions, where now the set L,t,s are deformation parame- 
ters and higher times. 

A fermionic representation for the tau functions will be written down and compared with the rep- 
resentation of appropriately deformed classical Dyson- Wigner ensembles, where an additional set of 
parameters s is added. The complex Ginibre ensemble will be deformed with the help of four sets of 
parameters, t,t',s, s', and related to the two-component Toda lattice. 



2 Gauss Real Ginibre ensemble and its deformation 

Real Hermitian ensemble. First let us remind that the ensemble of real Hermitian matrices (also 
known as Dyson- Wigner orthogonal ensemble) is given by the following partition function 

I d N -° E (t) = J d/i(X,t) , dn(X,t) = e-*tr(*')+trV(*,t) -Q dX .. (1) 

i>j 

oo 

V(x,t):=Y,x n tn (2) 

n=l 

where X is N by N real Hermitian matrix and and t = (t\, t2, ■ ■ ■) is a set of parameters which describe 
the deviation of the probability measure from the Gauss one. 

Let us restrict the space of our symmetric matrices to invertable symmetric matrices. Then more 
general deformation of the measure may be considered as follows 

i^- 0£ (I,t,s) = J dn(X,L,t t a), dn(X,L 7 t,s) = dctX L e- tTV{x ~ 1 ' s) dn(X,t) (3) 

where s and an integer L is the collection of new deformation parameters. Dyson- Wigner Gauss ensemble 
will be referred as G-OE, ensembles JTJ and ([3J respectively as d-OE and dd-OE. 

It was found by M.Adler and P. van Moerbeke in [7] that the ensemble ([1]) may be related to the 
so-called Pfaff Toda lattice [5]- Later in [5] J. van de Leur found that the orthogonal ensemble JT]) may 
be also identified with a tau function of the "large" BKP tau function introduced in [4] . Below we shall 
identify the ensemble © with a "large" 2-BKP tau function introduced in [3J. 

In terms of eigenvalues x\, . . . ,xm of the Hermitian matrix X the integral ([3J may be written as 

r N 
I d N d -° E (L,t,s)=a N / A( Xl ,...,x N )]lx?e-^+^)-vi*r\s) dxi (4) 

with some constant related to the volume of the orthogonal group O(N). Deformation parameters t 
and s are considered to be chosen in a way that the integral (HJ) is convergent. 

Real non-Hermitian ensemble. Let us turn to the so-called Ginibre ensembles. Gauss real Ginibre 
ensemble, also known as Gauss Ginibre orhtogonal ensemble (G-GinOE) is defined on the space of real 
matrices by assigning Gauss probability measure to each entry with the same variance: 

I%- G ™° E = J d»(X) , dn{X) = n e-WdXy = e-^ xx ") \{ dX tJ (5) 



2 



where X is N by N matrix with real entries. Measure (0 is invariant under the orthogonal transforma- 
tions of matrices X. 

The so-called elliptic deformation of G-GinOE measure as 

dn{X,a) = e- atlX ' 2 d^{X) (6) 

where a is the deformation parameter was studied in 1 . 

We shall consider the following deformation of the G-GinOE ensemble which we will refer as d-GinOE 
and which includes the elliptic deformation as a particular case 

Ip GinOE (L,t)= f dfi(X,L,t), dn(X,L,t) = dctX L e^^'V dfi(X) (7) 

where V is given by ([2J and where L = 0,1,2... and sets of numbers t = (t\, ta, . . .) are deformation 
parameters. We assume that t are chosen in a way that the partition function iff E (L, t) is finite. 

Remark 1. However even in the case where the integral is divergent itself, logarithmic derivatives with 
respect to the deformation parameters may be finite. How it works in the models of normal matrices, see in [BJ. 

Now let us consider the restricted Ginibre ensemble consisting of real invertable matrices and deform 
its measure as follows 

I d N d - GmOE {L,t,s) = J dfj,(X,L,t,s), dfi{X,L,t,s) = deUC L e tlV ^ x ' t) - tlV ^ x ~ 1 ^dn{X) (8) 

where X are real invertable N by N matrices and where V is given by ([5J and where an integer L and 
sets of numbers t = (t\, t2, ■ ■ •), s = (si, S2, . . .) are deformation parameters. This ensemble will be called 
double deformed restricted Ginibre orthogonal ensemble (dd-GinOE). Again we suppose that deformation 
parameters are chosen in a way which provides the existence of the integral lff'~ GmOE (L, t, s). Consider 
the following example. If one takes t — > t — Yli=i a i[Pi] an d s — > s — Yli=i namely 

tn -> tn OiP^ 1 , S n -> S n b 3l? ( 9 ) 

i=l j 

then the integral ((HJ reads as 

/ Y[det(l-p i X) ai J[det(l-q i X- 1 )~ h dv{X,L,t,s) (10) 
J * j 

Let prove that the integrals (0 and ([8]) may be identified with tau functions of the "large" BKP 
hierarchy and 2-BKP respectively. Then the size of matrices N and the set L,t,s play the role of the 
so-called higher times of this integrable hierarchy. 

To do it we should re- write integral (j5]) as the integral over eigenvalues. Let us remind that for the 
Gauss real Ginibre ensemble (namely, for the case t = s = 0, L — 0) this problem was solved in [2], see 
also a review article pQ. In very short it was done as follows. After applying the Schur decomposition 
X = U(A + AjW , where U is orthogonal, A is block-diagonal and A has nonzero blocks only above A 
they arrive at 

d^X) = e- 2Tl '( AAt+AAt >|£A||£A||n' (U-'dU)^ 

with the dashed product running over non-zero entries in the lower triangle of U'dU. The set of the 
eigenvalues Aj, i = 1, . . . , N of a real matrix X consists of real numbers and of complex conjugated pairs, 
denoted below as x% and (zi,Zi) respectively. Thus, the matrix A = diag(Ai, . . . , Ajv) may be written 
as diag(zi, z\, . . . , Zk, Zk, xi, ■ ■ ■ , %N-2k) with certain k. After integration over matrices A and U they 
come to the integral over eigenvalues only which includes integral over real values Xi and over complex 



2 It is also called the elliptic deformation of G-GinOE because in the large N limit the parameter a describes a 
deformation of a circle equilibrium domain of the eigenvalues in complex plane to an elliptic form. 
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eigenvalues Zj , Zj . Finely, after certain computations [5] the Gauss real Ginibre ensemble may be written 



as 



jG-GinOE = ^ . / A 2n (z ll Z 1 ,...,Z k ,Z k ,X 1 ,...,X N -2k)dng<m%_ 2k (11) 

k— "^^2fc,iV-2fc 

where the integration domain M.2 m ,N~2m is as follows Sftzi > • • • > 5ftz m , X2 m +i > ■ ■ ■ > xn, Zi e C+ 
(upper half-plane), a;^ € R, and where 

d0g 7l (z,z) = nerfc(^ 7 3 ) e^-V* , dfi^_ 2m (x) = [J e"**?** (12) 



with crfc(z) = -7= e 2:2 . Here the factor cn absorbs integrals over A and over U. Let us note that 
the Gauss case is a special case of ©: iG-GmOE = idf-GinOE( Q0Q y The factor ^ ig independent of 
L, t, s and was evaluated for GinOE in [2J. 

Now we notice that the deformation @ results in the multiplication of the measure dfi(X) by a factor 
which depends only on eigenvalues of the matrix X 

N 

dn{X) -> d/ipr, t, l, s) = dn(X) Y[ \L e v(\„t)-v(\-\ S ) (13) 

z=l 

where Ai are eigenvalues of X and V is given by (J5J). Then we introduce 

dng n -> dflgjt.L.s) = fjerfc (^^) l^l^e 2 *^.*)-^^- 1 ^-^^^ (14) 

JV 

dfi«_ 2m ^d^_ 2m (t,L, S )= J] xfey^.*)-^^ 1 .-)-^?^ (15) 

i=2m+l 

Let us introduce 
[f] 

Jn(L, t, s; a) = c* & / |z 4 | 2L A 2 „(zi, zi, . . . , z k , z k , xi, . . . , x N - 2 k) d£l C h {t, L, s) dQ§_ 2k (t, L, s) 

(16) 

We note that for a — > (up to a constant independent of the variables L, t, s) this expression is equal to 
the partition function of the deformed Dyson- Wigner orthogonal ensemble ([7]). When a — 1 this is the 
partition function for the deformed real Ginibre ensemble (|33l) . 

J N (L,t,s = 0;a = l) = b N I d N GmOE (L,t) (17) 

Jjv(L,t,s;a = 0) = a A r/^-° iJ (L,t,s), J N (L, t, s; a = 1) = b N I**- GinOE {L, t, s) (18) 

Independent of deformation parameters factors ajv and b^ may be found respectively in [10] and pQ. 

Tau functions and matrix ensembles. Let us remind the fermionic construction of the so-called 
tau functions [S] . Namely here we present a general fermionic expression for the 2-BKP tau function and 
relate it to ensembles ^ and (JSJ). In the next section it will be also related to a deformed symplectic 
Dyson- Wigner ensemble and to a deformed quaternion-real Ginibre ensemble. 

The very notion of tau function and it's fermionic construction was introduces by Kyoto school. Here 
we use it in a version suggested in j3] where an additional fermioc mode <j) was added, see Appendix lA.2l 

Following [5] we consider 

oo oo 

r + (t)=exp^i„^V,V4 + , i , r_(s)=exp^s„£>i</4-„ (19) 

n=l ieZ n=l igZ 

We shall need the following equality 

r+(t)^(z)r_(s) = C (t,s)e 1/ (^ t )- t/ (^ 1 ' s ))r_(s)^(z)r + (t) (20) 
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where 

oo 

c(t, s) = exp nt n s n (21) 

n=l 

The 2-BKP tau function may be presented in form of the following fermionic expectation value 

Ttr(L,t,a) = (JV + L|r+(t)e*r_(s)|£) (22) 

where $ is any quadratic expression in fermionic modes ipi,ip\, <p, see Appendix IA.2I Variables N, L, t, s 
plays the role of the 2-BKP higher times. Tau function (|22D solves 2-BKP Hirota equations, see [3] and 
Appendix IA.6I 

Remark 2 . The 2-BKP tau function is also a certain BKP tau function j4] with respect to the variables N, L, t 
and a certain BKP tau function with respect to the variables N, L, s (this explains the name 2-BKP). 

Tau function (|22p contains as a particular cases: tau functions of the Toda lattice hierarchy [S], [T5], the 
charged BKP tau function [4] and 2-DKP ("Pfaff DKP") tau function |20| . Each of these cases may be obtained 
by specifying the tau function (|22[) . namely, by specifying $ and s, see Appendix IA. 41 



Theorem 1. We have 

{-) NL c(t,s)J N {L,t,s;a) = (N + L\ T+(t) e a$c e* r r_(s) \L) (23) 

where 



z - z 



$ c = erfc [^pjr) ^(z)^{z)e-^ z d 2 z (24) 
is the integral over the upper half of complex plane and where 

$ r = ^J J sgn(xx - x 2 )il>(x 1 )'ip(x2)e~i x2l ~i x 2dx 1 dx2 + V2 J i/;(x)(j)e~i x2 dx (25) 

\NL jd-GinOEi 



Therefore thanks to (J7|) and Remark\M(-) N 1^ ( L ^) is a tau function of the "large" BKP while 

I d N d -° E (L,t,s) and(-) NL c(t,s)I d N d - 



thanks to (EH) both {~) NL c(t, s)I dd -° E (L, t, s) and (-) NL c{t, s)lff- GinOE (L, t, s) of the "large" 2-BKP 



ones. 



The sketch of the proof of ((23]) . We transform the vacuum expectation value in the right hand side to 
get the integral in the left hand side as follows. First we use Taylor expansion of the exponentials e Q * c 
and e^* r where only terms of the order a N ~ 2h /3 k , k = 0, . . . , [tH are non- vanishing between (N + L\ 
and \L). The term e Q * r should be considered in a way similar to [5] (see Appendix). Notice that the 
Taylor expansion of the exponentials of integrals yields the integration domain y^2k.N-2k when we re- 
write product of pairwise integrals as a iV-fold integral. Next we send T to the right and T* to the left 
taking into account relations (N + L\T*(s) = (N + L\ and T(t)\L) = \L) and using ([2U)l . In this way 
we get rid of operators Y and Y* and instead obtain factors e v inside integrals, the integrals are still 
between (N + L\ and \L). At last we get rid of fermions thanks to 

(N + L\tp(zi) ■ ■ ■ ip(z N )\L) = Zi ■■■ z L tA n (z) 1 A n (z):= \{{z l -z j ) (26) 

i<j<N 



obtained by a simple direct calculation. We obtain 

Now the fermionic expectation value (TV + L\ Y(t) e Q * c e* r r*(s) \L) is an example of tau function 
with respect to the variables N,L,t introduced in j4] where Hirota equations for such tau functions were 
written down. By symmetry it is also tau function with respect to the variables N, L, s. Such ("coupled" 
BKP, or, the same 2-BKP) tau functions were considered in [3]. 

Thus we establish that the partition function of the deformed quaternion-real Ginibre ensemble ([7]) 
is the subject to the theory of integrable systems. 
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3 Quaternion-real Ginibre ensemble and its deformation 

Each TV x TV quaternion matrix (or, the same, a matrix with quaternion entries) may be viewed as 
27V x 27V matrix if the quaternions e„, n = 1, 2, 3, 4, are realized as 2 x 2 matrices: 

e o = I I ? J . ei = ( * _°, J , e 2 = [ ° J. J , e 3 = ( ° * j (27) 



Then quaternion-real matrix is a quaternion matrix where each entry is a linear combination of quater- 
nions with real coefficients. Being viewed as 27V x 27V matrix a quaternion-real matrix X may be 
recognized via the property 

EXE = -X (28) 

where the bar means complex conjugate and where E is the block-diagonal 2 TV x 2 N matrix with matrices 

^ ^ on the main diagonal; this follows from the explicit expression (|27l) and from e2&j&2 — Sj, 

j = 1,3 and from e2ej&2 = —6j, j = 0,2. Below we shall treat quaternion-real matrices as 2N x 27V 
matrices with the defining property (|28[). 

Quaternion-real Hermitian ensemble. First we recall the case where in the quaternion-real 
matrix X is Hermitian (in this case X is called self-dual). The ensemble of random self-dual matrices 
is called symplectic ensemble, the Gauss symplectic ensemble (G-SE) is called Dyson- Wigner symplcctic 
ensemble. We will consider a deformed symplectic ensemble (d-SE) as follows 

I d N SE {L, t) = J dn(X, L, t) , dfi(X, L, t) = e -* 2 +^v(x,t) d/i pQ (2g) 

where dfi(X) is the measure on the space of self-dual matrices (for details see [TO])- As in the orthog- 
onal ensemble case we may consider the subspace of intertable matrices and consider double deformed 
restricted ensemble of quaternion-real matrices (dd-SE) 

I^ d - SE (L, t, s) = J dn(X, L, t, s) , dn(X, L, t, s) = d ctX L e- a;2+2trV '^ t ^ 2trV '( x_1 ^ dfi{X) (30) 

where we add the deformation parameters s = (s2, s%, . . . and integer L. Being re- written as an integral 
over eigenvalues of Hermitian quaternion real matrices X |10] it is as follows 

lff~ SE (L,t, S )=d N j {^{ Xl ,x 2 ,...,x N )f\{x} L e-^ v ^-^ 1 ^U^ (31) 
•* »=i 

where jci, x\, X2, xq, . . . , xn, %n are eigenvalues of the self-dual random 27V by 27V matrix X and the 
factor djv does not depend on deformation patameters. 

Quaternion-real non-Hermitian ensemble. Then Gauss quaternion-real Ginibre ensemble, also 
known as (Gauss) Ginibre symplectic ensemble (G-GinSE) is defined on the space of quaternion-real 
matrices by assigning the same Gauss probability measure to each entry: 

jG—GinSE = J = e -itr(^t) J| ^. . ( g 2) 

where X is treated as 27V by 27V matrix. Measure (13"2l is invariant under the symplectic transformations 
of matrices X. 

Remark 3. A Hermitian quaternion-real matrix is called quaternion self-dual matrix. Ensemble (|32p is non- 
Hermitian analogue of the ensemble of random quaternion self-dual matrices (known as Dyson- Wigner symplectic 
ensemble), see [10] for details. 

The elliptic deformation of G-GinSE measure [I] is quite similar to ©: d^(X) -> e ^ atr ( x +<- x ) d^(X) 
where a is a deformation parameter. 

The partition function for the deformed quaternionic-real Ginibre ensemble (d-GinSE) will be defined 

as 

Ip GinSB {L, t) = J dfx(X, L, t), dfi(X, L, t) := detX L d^{X) (33) 



G 



where a non-negative integer L and the set of numbers t = {t\,t2, ■ ■ •) are deformation parameters. 
On the space of invertable quaternion-real matrices we consider double deformed Ginibre ensemble 

lf- GinSE {L,t,&)= J dn(X,L,t,B), d^(X,L,t,s)=detX L e iTV ^- tlV ^ x ~ 1 ^d^(X) (34) 

where an integer L and and the sets t = (ti, t 2 , ■ ■ •), s = (s\, s%, ■ ■ •) are deformation parameters. Here 
we assume that the integral is either convergent or regularized, see Remark [1] 
Notice that if one takes 

t n -> t n - — a 'P"' s ™ ^* s " ~ 2n X! b i q j ( 35 ) 

i=l ~ j 

then the integral (l33|) reads as 

/ detX L Y[dct(l-p l X) a ''[[dct(l-q l X- i y b ' dfi{X,L, t,s) (36) 

As in the case of real matrices all complex eigenvalues occur in complex conjugated pairs. Real 
eigenvalues have multiplicities no less than two. A calculation similar to calculation for real matrices 
yields [T] (see also [27], [21] , QH]) 

lff- GinSE (L,t,s)= (37) 
A 2n (z u z u Z2, z 2 , . . . , z n , z n ) \zi ~ zi\ | Zi |M e v^,t)+v(s ( ,t)-v(, r \.)-v( Sr S.)-| ft |» dazj 

with some constant c^r where the integration domain M n consists of the sets of z where $izi > $lzi + i 
and Qzi > 0, and V is given by $Z§. 

Theorem 2. Introduce 

c(t,s)Jjv(i,t,s,a,/3) = (2A r + L|r + (t)e Q * c e /3<E> ''r_(s)|A0 (38) 

w/iere 

/ (z- z)V'(^)V'(^)e~ |z|2 d 2 ^ (39) 



<1>; = / ^^^{x)e- x ' dx (40) 



is i/ie integral over the upper half of complex plane. Then 



I d - GinSE {L, t) = d N J N {L, t, s = 0, a = 1,0 = 0) (41) 



lff- GmSE (L, t, s) = aV Jat(L, t, s, a = 1, = 0) , I dd ~ SE (L, t, s) = c N J N (L, t, s, a = 0, (3 = 1) (42) 



Therefore Iff (L, t) is an example of the large ("fermionic") BKP tau function of J^j with respect 
to the variables N,L,t. Then c(t, s)lff' GinSE (L, t, s) and c(t, s)I% l ~ SE (L, t, s) are examples of the 
large 2-BKP tau function considered in with respect to the variables N,L,t,s and of large 2-DKP tau 
function introduced in [20], see Appendix \A.4\ 

Let us note that the presentation with s = 0, namely, the representation for the symplectic Dyson- 
Wigner ensemble (where a = 0) was found by J. van de Leur in [S]. 

3.1 On perturbation series in deformation parameters 

In [5] we have considered the following series over partitions 

J2 A h(x) (L)8 X (t) (43) 

Aep 

e(\)<N 
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where s\ is the Schur function (see Appendix IA.3|) and j?/i(a) [L) is a certain Pfaffian, see Appendix 
IA.5I Such series one obtains for many kinds of matrix integrals, see |29j . for instance for integrals over 
symplectic and orthogonal groups, see Appendix IA.7I The notations are as follows. A = (Ai, . . . , Xn), 
Ai > . . • Ajy > is a partition, see [11] and h(X) = Xi — i + N are the so-called shifted parts of A. The 
factors Ah on the right-hand side of (|43|) are determined in terms a pair (A, a) =: A where A is an infinite 
skew symmetric matrix and a an infinite vector. For a set h = (h±, . . . , /ijv), h± > • •• > the numbers 
Ah are defined as the Pfaffian of an antisymmetric 2n x 2n matrix A as follows: 

A h (L) := Pf[i] (44) 

where for N = 2n even 

Aij = -Aji := A hi+Lthj+L , l<i<j<2n (45) 

and for N = 2n — 1 odd 

hi = -A* '= I Ah ^ +L ., ^ I ~ 1 K j f 2n 1 • (46) 
3 3 \ a hl +L if l<Kj = 2n v ' 

In addition we set Aq = 1. 

Having the fermionic representation it is straightforward to write down the perturbation series of 
partition functions of Ginibre ensembles in the Schur functions like it was done in other cases (for 
instance see [23]. fig]. [17] ). 

(I) First let us do it for the quaternion real Ginibre ensemble. Re-writing 

f (Z - ^{z^e-^e-^^^^dH = AnmMm (47) 
"' c + n,m 

where A nm is a moment matrix 

A^ SE (s)= [ z n z m (z-z)e~^ 2 - v ^ 1 ^- y ^ 1 ' s U 2 z, af nSE = 0, n,m = l,...,N (48) 
Jc+ 

we reduce problem to the solved in [3]. We obtain the following series in the Schur functions (see [TT] 
for definitions and details) s\ (please, do not mix with the deformation parameters Sj): 

lff- GinSE (L,t, S )= Afhf E (L, s) Sx (t) (49) 

A 

1(\)<2N 

In case s = we get A m , n = -A n?m = n!J m +i jn . 

(II) For DW symplectic ensemble the formula is the same, but now instead of 1|48)) we have 

A s n E (s) = ^ [ x^-^-^^dx, a s n E = Q, n,m = l,...,N (50) 

(III) For the DW orthogonal ensemble we have in case N is even the formula is the same, However 
for N odd we need to define A^ h }- 

A™(s)= [ [ x n y m sga(x-y)e-i x2 -^ 2 - v ^ 1 ^- v ^~ 1 ^dxdy, (51) 

JR JR 

flf (s) = f x n e~ v ^ 1 ' s Ux (52) 
Jr 

(IV) For the real Ginibre ensemble we have 

A°%° E {s) = j z n z m eric (^jf) e-^-^'^-^^z (53) 

+ f [ x n y m sgn(x-y)e'^ x2 -^ 2 - v ^ 1 ^- v ^ 1 ^dxdy, 
Jr Jr 
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a 



GinOEl 



(s) = / x n e~ v ^ X > s) dx (54) 



For TV even (which occurs in the cases of DW orthogonal and in real Ginibre ensembles) the formula for 
Ar/j} is slightly more spacious, see [3]. 

Let us also note that in all cases where s = (namely, for d-OE, d-GinOE, d-SE, d-GinSE ensembles) 
moments can be explicitly evaluated in a straightforward way. 

4 On Pfaffian formulae 

Various Pfaffian formulae are known in the study of non-Hermitian matrices, see |26j . PQ. These formulae, 
and perhaps some new ones, may be obtained by applying the Wick's rule to the evaluation of the 
fermionic expectation values, see Appendix IA. 51 

Example. Let us choose t as t shifted by the variables © where all pi are different, all bi — 0, all 
Oi = —1 and Ni — N even. Then according to relation and the Wick's rule one can write the 
expectation value as a Pfaffian of the pairwise expectation values, see Appendix IA. 51 Then (|10|) is 



f N n AT (L + 1)(2-JV) 

J J] dct (1 - MX, L, t, s) = Uj^L——^ Pf [K] n>m=1 _ N 



(55) 



li>j 

where the fermionic representation for the tau function tn(L, t,s) (l2"2l yields 

K nm = (L|V' t (Pn)'0 t (Pm)e*|L) = {p m ~ Pn)K* nm (56) 

where K* m . 

K* nm = x L y L l - V]e n r-^y- forOE (57) 

Jr (1 - xp n )(l - yp n )(l - xp m )(l - yp m ) 



erfc 



( \ z ~ z \ \ e -3iz 2 + V(z,t) + V(z,t)-V(z- 1 ,s)-V(z- 1 ,s) ( j2 z 

Km= I (z~z)\z\ 2L — rjz — (58) 

Jc+ (1 - zp n )(l - zp n )(l - zp m )(l - zp m ) 

+ / x L y Ll -^-^ m m m forGinOE (59) 

Jr (1 - xp n ){\ - yp n )(l ~ xp m )(l - yp m ) 

r -x 2 +2V(x,t)-2V(x- 1 ,s) j 

K ™ = / x2L a X v Wi X v v forSE (60) 

K ™ = L '1 - (i-,p..)d -»..)d foi Gi ° sE (61) 

Here z is the complex conjugated to z. In this example the number of parameters pi was equal to N. In 
other cases we typically obtain pfaffians of block matrices. This will be written down in a more detailed 
text. 



5 Gauss complex Ginibre ensemble and its deformation 

The complex Ginibre ensemble is considered to be easiest and relatively studied and it is quite similar 
to the well-known model of normal matrices [16 , therefore I will skip details. 

Hermitian random matrices. First we remind that the ensemble of random Hcrmitian matrices X 
with Gauss distribution for each entry was is known as Gauss unitary ensemble (GUE). The deformation 
of the measure caused by multiplying the Gauss factor e~~ trX by e tTV ^ x,t ' was related to integrable 
systems in [T2]. Let us added additional deformation parameters L,s: 

J det(X) L e^^ 1 (( t ™ +t -) trXm -( s - + <") trf_m )- trX2 d / i(X) (62) 
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where d^i(X) is the Haar measure on the space of Hermitian matrices (see [TU] for details). Actually we 
do not need t' and s' but we keep it to compare the model with the non-Hermitian one. The fermionic 
representation for this model (with s = s' = 0) was given in [TJ5] (see also [H]). However to compare the 
result with the non-Hermitian case it is suitable to use another fermionic construction considered in |14) 
as follows 

e N (L + N, -N\ r(t, t') e Ic^ ll) (x)^ m (x)e-* 2 dx p (Sj B /j | Xj Q> (63) 

with certain independent of deformation parameters factor ejv- Here we add the additional deformation 
caused by L, s. 

Non-Hermitian random matrices Complex restricted Ginibre ensemble with the double deformed 
measure 

J det(X) L Met(Xt)-^ e ^=i( i ™ texm + t ™ tr ^ t ) m -^trx-"'-^tr(xt)--)-trxxt Y[ dX .. (64) 

(where parameters L\, t, t', s, s' are chosen in a way that the integral is either convergent or may be 
regularized (see Remark [T])) is equal to the following two-component 2D Toda lattice [9 tau function 

f N {Lt +N,L 2 - N\ T(t, e/c V- (1) (^ t<2) (*) e -wV* f ^ ^ ^ ^ (65) 

The factor /jv is independent of L\, L 2 , t, t', s, s'. 

This type of representation was previously used in |17) in the context of two-matrix and normal 
matrix models. Perturbation series in deformation parameters for the complex Ginibre ensemble are 
basically the same as found in |17) . 
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A Appendices 

A.l Vertex operators 

Vertex operators we need are as follows 

X(L,t,X) := e £~=i^ A £ e -£r=id^4r , x\L,t,X) := e~ s ~=i A "*™ X~ L e E ^ (66) 

Y(L,s,X) := e -^^ x ' ns "X L e^=^^ , Y\L,s,X) := e^"° =1 x ™ Sn X~ L e~ ^«=i (67) 

Interesting historical fact that the formula which relates fermions to bosons first was found in [22 . 
The following bosonization relation is useful 



ii + wir ft + Vr» n - Wfrr'i-'^'Wff+W „„ 



Introduce 

n n (L,t)=resf ^; t ' A) lt (jL)t ,A) N ) , h* n {L, s) = res (ft (L, s, A) ) (69) 



A. 2 Fermions 



We shall remind some facts and notations of [5]. Introduce free fermionic fields ip{z) = 'Yln^'tyiZ 1 , 
^/t (z) = X^iGZ ^-i-i z% wnose Fourie components anti-commute as follows ipiipj +ipj4'i = 4*1 Yj +V- , j'</'i = 
and tpiipj + i/jji/Ji — Sij where Sij is the Kronecker symbol. We put 

^|0) = ^|0) = (0|V-< = = (70) 

where (0| and |0) are left and right vacuum vectors of the fermionic Fock space, (0| • 1 • |0) = 1. Also 
introduce 

'''IW-!-^ if n<0' ln) -\^ n --^U\0) if ™<0 [n) 

Then (n| • 1 • |m) = <5„, m . 

Following [3] we introduce an additional Fermi mode which we shall denote by <fi with properties d 

<M+^0 = #J+^ = O, 2 -^ (72) 
#)} = |0)-^, W=i<0| (73) 

such that (L\(j)\L) = 

Now, two-component Fermi fields used in Section [5] are defined as 

^)(z) = ]>>> 2 „ +l , ^ wt (*) = E*"""VL + i ( 74 ) 



3 In notations of [4] our ip n , ipX and </> read respectively as i/> n ii, i an d V'O 
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Other details about two-component fermions may be found in [S] or in [17] . 
We have 

X(L,t,X)X^L,n)(N + L\T + (t)gT_( S )\L) = (N + L\T + (t)^(X)^ (^gT_(s)\L) (75) 
yHl,s,^Y(X)(N + L\T + (t)gT_(s)\L) = (N + L\T + (t)grP(X)^ (fi)T_(s)\L) (76) 
Then it follows that 

Q n (L,t){N + L\T + (t)gT_(s)\L) = (N + i|r+(t)fi n5 r_(s)|L) (77) 
Q n (L,t){N + L\T + (t)gT_(s)\L) = (N + L|r+(t)gf2 n r_(s)|L) (78) 

where 




Using the fermionic representation one may verify that tau functions related to the deformed orthogonal 
and deformed symplectic ensembles obey the constraints 

(o„(L,t)-n*(L,s))r(L,t / ,s / ) = 0, n>l odd (80) 

where t' k = t% — \&2,k s' k = s fe — \&2,k (this shift appears due to the Gauss measure in undeformed 
ensembles) . 

A. 3 The Schur function 

Consider polynomials h n (t) defined by e^™ =1 z *™ = Y^=o z ™^n(t). Then the Schur function labeled by 
a partition A = (Ai, . . . , A& > 0) may be defined as s\(t) = dot (h\ i -i + j(t)) i - =1 fe . 

A.4 From 2-BKP to TL, BKP, 2-DKP, DKP 

The general expression for $ of (l2"2"j) is as follows 

ijez j,jez ijez iez iez 

To get TL tau function [5], [T^] we put all Aij,Bij,ai, hi and AT to be zero. 
To get BKP g] we put s = 0. 
To get 2-DKP [20] we put all a, = b t = 0. 
To get DKP [9] we put s = and all a, = b, = 0. 

Tau functions (f2"3"]l and (j3"8"]) correspond to the case where all Bij , Dy , &j vanish (and for ((55]) also all 
aj = 0). The case where only Ajj (and perhaps a,) are non- vanishing is characterized by the condition 
that the so-called wave functions 

w^(N,L,t,s,X) = *(WW(M,s) = x L e v ^- v ^~ l ^P N (L,t,B„X) (82) 

T N [L,t,s) 

w (°)(7V,L,t,s,A) = rf(L ' S ' A r )TJv( f' t ' s) = A £ e v ( A > t )-^ A - 1 '")Q w (i > t,8 ) A) (83) 

t n\L, t, s) 

are equal, and Pjv = Qn are polynomials in A of the order N. 
A. 5 Pfaffians 

If A an anti-symmetric matrix of an odd order its determinant vanishes. For even order, say k, the 
following multilinear form in Aij ,i < j < k 

Pi i A ] := X! S S n ( Cr ) 4t(1), CT (2) At(3),<t(4) ' ' ' K{k-l),a{k) ( 84 ) 

(7 

where sum runs over all permutation restricted by 

a : a(2i- 1) < a(2i), a(l) < <r(3) < ■ ■ ■ < a(k - 1), (85) 

coincides with the square root of detA and is called the Pfaffian of A, see, for instance [TUJ- As one can 
see the Pfaffian contains 1-3-5 (k — 1) =: (k — 1)!! terms. 



13 



Wick's relations. Let each of Wi be a linear combination of Fermi operators: 

= ^ V im1pm + U im^ln > * = 1, • • • , " 



Then the Wick formula is 

(l\w 1 ...w n \l) = { P Jl A kj=i,...,n ifniseven 
1 otherwise 

where A is n by n antisymmetric matrix with entries Aij — (l\ibiWj\l) , i < j. 



(86) 



A. 6 Hirota equations 

Hirota equations for the large BKP hierarchy were written in [J] . For 2-BKP hierarchy Hirota equations 
are as follows [5] 



dz 
2nd 



„N+l-N' -I' -2 „V(t-t' iii u , r_-l. ii + r,-li c 

-— z e ■ >T N i + i{l ,t + [z \,s )TN-i[l,t— [z ,s 

2m 

= I ^/-^ v ^'- s ^t n ^ 1 (1' + l,t',s' - [z])t n+1 (1 - l.t.s- [z] 
J 2m 

+ I ^-V^"Vi(l'-l,tV+[«]) W .i(i + l 1 t,S+[2] 
J 2m 

(-!)''+« 



(l-(-l) iv +iv )r^(r,t',s')r Ar (/,t,s 
The difference Hirota equation may be obtained from the previous one [3] 







Tjv(i, t + \p- x ])T N+1 {l, t + [a" 1 ]) - -2— TN (l, t + [a-^Jrv+i {I, t + [/3" 1 ]) 

p — a 



a- (3 



+ ^T N+2 (l,t + [a- 1 ] + [/T^rjv-i^t) = Tjv+i(/,t + [a" 1 ] + [/T^Tjv^t). 

A. 7 Integrals over orthogonal and symplectic groups 

Using 



/ s x (0)d*0 = \ 

JO£<0{N) 1 



1 A is even 
otherwise 



seSp(N) 



s x (S)d*S={ 1 



X tr is even 
otherwise 



we get 



N):= ( e ^^° m d*0= V SA (t) 
JoeO(Af) 



J 2 (t,N) 



Jsen P (2n) lt , ______ 



«(A)<2n 



(87) 



(89) 
(90) 
(91) 



The right hand sides were obtained in [5] as examples of the BKP tau function. Thus integrals 
Ji(t,N) and J2(t,N) are tau functions, N and t being BKP higher times. 
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